This work had the objective of identify existing climate data and coffee leaf rust incidence data for the state of Oaxaca, Mexico. The initial approach was to gather as much information available, and evaluate if it was enough to identify correlations between leaf rust incidence and climate variability. All the information gathered has been stored in an online Github repository, along with the R and MatLab code used to download, clean, and join all the information. Additionally, ESRI shapefiles for Mexico’s states and municipalities are provided. The repository lives permanently at https://github.com/jcvdav/OaxacaLeafRust. Information can be downloaded by downloading a *.zip folder. The entire folder weights about 0.5 GB.
Databases are provided as Excel spreadsheets (*.xlsx), MatLab files (*.mat), R data (*.RData) and comma-separated values (*.csv). The Excel files contain a spreadsheet with Metadata, describing each field (column) of the database, as well as the source of the information and relevant comments. Climate data from CLICOM (Deliverable 3) is not available in Excel format due to its large size.
Based on the observed spatial distribution of data for Oaxaca, it is not possible to pursue further correlation analysis between climate data and coffee leaf rust incidence. The closest climate station-leaf rust data pair is that of the municipality “Candelaria Loxicha”, with a station in the neighboring municipality of “Pluma Hidalgo”. However, they are approximately 50 km away along a mountain range, which will likely causes environmental conditions to be very different. Instead, exploring the neighboring state of Chiapas might yield some more information.
A literature review on the most influential papers on coffee leaf rust and environmental variables was performed. From this effort, the most relevant climate and management factores are included in Table I. It is important to mention that the identification of management strategies is critical, as the direction and magnitude of interactions between climate variables and coffee leaf rust might be affected.
Bock, K.R. 1962. Seasonal Periodicity of Coffee Leaf Rust and Factors Affecting the Severity of Outbreaks in Kenya Colony. Trans. Brit. mycol. Soc, 45(3): 289-300.
Avelino, J., Willocquet, L., Savary, S. 2004. Effects of crop management patterns on coffee rust epidemics. Plant Pathology, 535: 541 – 547.
Avelino, J. Zelaya, H., Merlo, A., Pineda, A., Ordoñez, M., Savary, S. 2006. The intensity of a coffee rust epidemic is dependent on production situations. Ecological modeling, 197: 431-447.
Avelino, J., Romero-Guardián, A., Cruz-Cuellar, H.F., Declereck, A.J.F. 2012. Landscape context and scale differentially impact coffee leaf rust, coffee berry borer, and coffee root-knot nematodes. Ecological Applications, 22(2): 584 – 596.
Avelino, J., Cristancho, M., Georgiou, S., Imbach, P., Aguilar, L., Bornemann, G., Läderach, P., Anzueto, F., Hruska, A.J., Morales, C. 2015. The coffee rust crises in Colombia and Central America (2008 – 2013): impacts, plausible causes and proposed solution. Food Security, 7: 303 – 321.
-Cilas. C., Goebel, F.R., Babin, R., and Avelino, J. 2016. Chapter 6: Tropical Crop Pests and Diseases in a Climate Change Setting – A Few Examples. Torquebiau (ed), Climate Change and Agriculture Worldwide. DOI: 10.1007/978-94-017-7462-8_6
This contains daily data for climatic variables collected from weather stations at 2 m above the ground. Data comes from INIFAP’s website at http://clima.inifap.gob.mx/redinifap/, where one can search by state, station, year, and month. This same portal has station data available for all other states, including Chiapas, Veracruz, and Guerrero, states that also produce coffee. Data were obtained by copy-pasting directly from the tables displayed in each query for the state of Oaxaca and every station listed under it. There are a total of 49 stations with data that covers some point between 2007-02-23 and 2016-10-19 (Fig. 1). Within that window, some stations may have gaps of missing data. The compiled database can be downloaded from github in *.csv format or in *.RData format. This data has the following fields:
Name: the name of the station
Municipality: the municipality where the station is
Latitud1: Latitude in decimal degrees
Longitude1: Longitude in decimal degrees
Latitude2: Latitude in DMS format
Longitude2: Latitude in DMS format
Date: Date in dd//mm/yy format (*.csv) and YYYY/mm/dd (*.Rdata)
Preci: Precipitation in mm
Tmax: Maximum daily temperature in °C
Tmin: Minimum daily temperature in °C
Tmed: Mean daily temperature in °C
VVMax: Maximum wind speed in km/h
DVVMax: Direction of the maximum wind speed in azimuth
VV: Mean wind speed in km/h
DV: Mean direction of wind speed in azimuth
RadG: Global radiation in w / m^2
HR: Relative humidity (%)
ET: Reference evapotranspiration (mm)
EP: Potential evapotranspiration (mm)
Coffee leaf rust data weekly reports from royacafe. The parent directory for all the reports is available at http://royacafe.lanref.org.mx/ReportesSPEyC_doc/. The weekly reports include information for coffee leaf rust incidence. The data is reported at municipality level, and no specific coordinates or information about the number of coffee plots is provided by the source. This data set contains information for a total of 13 municipalities, with information from 2015-05-04 to 2016-08-21 (Fig. 2). There are three main clusters of coffee leaf rust information. One of them is located in the Northernmost part of the state. For that region, there are between 1 and 5 weekly reports for each of the 4 municipalities. The second cluster is in the Eastern region of the state, where 3 municipalities have 1 observation each. The third cluster is located in the South of the state. This cluster has 6 municipalities, with between 6 and 28 observations.
This data set contains information for:
Municipality: Municipality
Initial: Initial date of sampling for the weekly report as day/month/year
Final: Final date of sampling for the weekly report as day/month/year
PS: Plant severity, expressed as percentage of affected plants
PVI: Plant variability index, expressed as relative variability from regional mean
LS: Leaf severity index, expressed as percentage of affected leaves
LVI: Leaf variability index, expressed as relative variability from regional mean
AL: Number of affected leafs per plant
Lechoso: Number of stage-1 fruits per sampling site
Consistente: Number of stage-2 fruits per sampling site
Maduro: Number of stage-3 (mature) fruits per sampling site
Note: While these are called weekly reports, there are not 52 reports per year. Apparently, this makes reference to the fact that the sampling period for each report is often a week.
This data comes from CLICOM’s database at clicom-mex.cicese.mx. The portal allows the user to download sets of 20 files (stations) at a time, either in *.csv or *.mat formats. This contains daily data for climatic variables. There are a total of 162 stations with data that covers some point between 2000 and 2012. The oldest station has data from as far back as 1922. While this data has a better spatial distribution (Fig. 5), information for recent years is not available.
This data set contains information for:
id: Station identification number
station: Station name as “STATION-STATE”
latitude: Latitude in decimal degrees
longitude: Longitude in decimal degrees
day: Day
month: Month
year: Year
date: date as string of characters in day-month-year format
datenum: MatLab’s datenum
id2: Combination of station id and datenum as unique identifier for every day
Prec: Precipitation in mm
Evap: Evaporation in mm
MeanT: Mean daily temperature in °C
MinT: Minimum daily temperature in °C
MaxT: Maximum daily temperature in °C
All figures are available as interactive figures at
Figure 1 - Initial and final times of sampling for each station. These window of meassurements may be interrupted at different points for different stations.
Figure 2 - Initial and final times of coffee leaf rust data. These window of meassurements may be interrupted at different points for different stations.
Figure 3 - Map of Oaxaca and municipalities. Color of municipalities represent the number of weeks (weekly reports) with available information for each municipality, also shown as text over each municipality. Red dots indicate the locations of weather stations from which data is available from INIFAP’s database.
Figure 4 - Initial and final times of sampling for each station. This represents a subser of all stations to reduce the size of the graph. This subset represents stations with more than 30 years of data, from 2010 back. These window of meassurements may be interrupted at different points for different stations.
Figure 5 - Map of Oaxaca and municipalities. Color of municipalities represent the number of weeks (weekly reports) with available information for each municipality, also shown as text over each municipality. Red dots indicate the locations of weather stations from which data is available up to 2012 from CLICOM’s database.